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We present an approach for computing the real parts of the nonleptonic B — >■ DP and B DP 
{P = K, tt) decay amplitudes by using lattice QCD methods. While it remains very challenging to 
calculate the imaginary parts of these matrix elements on the lattice, we stress that their real parts 
play a significant role in extracting the angle 7 in the b—d unitarity triangle of the CKM matrix. 
The real part on its own gives a lower bound to the absolute magnitude of the amplitude which is in 
itself an important constraint for determining 7. Also the relevant phase can be obtained by using 
B-decays in conjunction with relevant charm decay data. Direct four-point function calculations 
on the lattice, while computationally demanding, do yield the real part as that is not impeded 
by the Maiani- Testa theorem. As an approximation, we argue that the chiral expansion of these 
decays is valid in a framework similar to that of hard-pion chiral perturbation theory. In addition 
to constructing the leading-order operators, we also discuss the features of the next-to-leading order 
chiral expansion. These include the contributions from the resonance states, as well as the generic 
forms of the chiral logarithms. 

PACS numbers: 11.15.Ha,11.30.Rd,12.38.Gc,12.39.Fe,12.39.Hg 



I. MOTIVATION AND INTRODUCTION 

The B+ D°P+, B+ ^ D P+ decay processes (with P ^ K,tt) are of significant phenomenological importance. 
These decays can be used for a direct, data-driven, extraction of the CP-odd phase 7 of the h—d unitarity triangle 
in the CKM matrix. In principle, given sufficient number of _B-mesons they can provide a determination of 7 from 
experiment to an unprecedented accuracy of « 0.1% [l|. In comparison, the projected accuracy of the angle /3 is at best 
« 0.5%, and for a it is likely limited to a few percent [^l- Using only the charged B meson decays to achieve a precise 
determination of 7 is highly valued, since the underlying decay modes are dominated by tree-level weak-interaction 
processes. 

These direct methods for deducing 7 from charged B meson decays involve interference between and decays 
to common final states,^ for example 

B- D°P- fp- and ^ ifp^ fP^ , (1) 

where / = i^gTr^, JC+if ~, isr~7r+, A'~p+, A'*~7r+, etc. 0-0I- particular, when the final state / is such that decays 
from D are Cabibbo allowed, but those from are doubly-Cabibbo suppressed, e.g. A'+tt", then in the overall 



charged B decays, there is a tendency for the interference to be maximal, giving rise to the possibility of large 0(1) 
CP-asymmetries [5|. This is important as larger asymmetries tend to require fewer B mesons for deducing 7. 
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^ Note that this involves — 'D*' mixing and in the Standard Model, CP violation from this source is assumed to be very small; this is 
assumed in all 7 analyses so far [J]. 
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Since the time these methods were proposed^it has always been recognized and emphasized that studies of charm 
decays can be very helpful for extracting 7 p, H, H, in particular, precise knowledge of the branching ratio of 
the relevant charm decay modes and the strong phase(s) can significantly facilitate the determination of 7. Specific 
methods [lo| have been proposed for studies at charm facilities for this purpose and great deal of experimental activity 
has taken place and progress is being made [ll| - [l5| . 

The methods proposed in Refs. 0, Q allow the extraction of 7, as well as the relevant strong phase difference in 
D^P^ and B^ D amplitudes. In these methods, the branching ratio, 

Br[B- D"P-], (2) 

is an essential input. This branching ratio (for P — K or tt) has been experimentally measured with good precision. 

Due to technical reasons, Br[i?^ — > D^P^] is not accessible to experiment Q; for this reason in the method of 
this branching ratio, expressed as the ratio, 

= Br[i?- ^ D^P-] ' 

is treated as an unknown that can be solved for along with 7. However, determination of this ratio in addition to 7, 
places additional demands on the number of B mesons that are needed. For this reason, despite the large statistics 
of the two S-factories 0(10^) B meson samples], 7 is presently determined to only ~ 0(25%). This should be 
compared with about 3% for /3, and about 5% for a. To further improve the accuracy on 7, inputs from lattice QCD 
(LQCD) on the ratio in Eq. ^ would be very useful. In other words if the lattice could provide an accurate value of 
this ratio, a fewer number of B meson samples will be needed to achieve a given accuracy on 7. 

For the purpose of a lattice study, we define a "reduced" ratio which is independent of the CKM matrix elements, 

red = rBP _ IVelKgP , . 

'BP — T/combo ~ IT/* 1/ 12 ' 

^CKM l^ub^cql 

where q ~ s,d depending on whether P — K or t:. Needless to say, the s tudy of hadronic weak decays on the lattice 
continues to represent an outstanding challenge. Exploratory studies [16l - [l8| initiated in the 80's did not have much 
success, because of the Maiani- Testa no-go theorem (MTNGT) [l^. This theorem states that Euclidean four-point 
correlators (three sources for external hadrons plus one weak-operator insertion point) always result in the average 
of in- and out-states, leading to the impossibility of extracting information about the strong phases. That is, one 
can only compute the real parts of nonleptonic decay amplitudes from such correlators in a finite volume.^ For the 



calculation of if — >■ tttt on the lattice, one can evade the MTNGT using the Lellouch-Liischer (LL) method [2l|, and 
the RBC-UKQCD collaboration is making considerable progress [22|, [23| in this avenue. 

On the other hand, the lattice computation of these B DP decay amplitudes remains challenging, both because the 
calculation of the lattice four-point function is computationally demanding to evaluate and because the LL method is 
only applicable to processes involving elastic final-state scatterings. With the advent of new powerful computers such 
as the BG/Q, the former difficulty may be overcome in the near future, especially since lattice results for the real 
part of these amplitudes could provide valuable information on the ratio Eq. and thereby help in the extraction 
of 7 when combined with experimental measurements. For now, we will investigate the use of some approximation 
methods for tackling these amplitudes.'^ 

Let us recapitulate that for a determination of the ratio, Eq. Q, what one needs is the absolute value of the amplitudes 

for the B — )- D'^P and B D P modes, and not just the real part of the amplitudes that is accessible on the lattice. 
So what we envision is that the phase of these amplitudes will also be accessible by combining information from the 
method of js], Q with the phase of the relevant charm modes coming from charm studies as briefly alluded to above. 



^ The D—P spectrum in finite volume is rendered discrete, enabUng the extraction of the energy of the excited state which corresponds 
to the physical state [20II . 

^ It is useful to note that lattice calculation of these B — > D{D)P amplitudes involve no mixing with lower dimensional operators, 
"eye-graphs" or disconnected diagrams and to that extent are simpler than K ^ 2Tr amplitudes in the A7 = 1/2 channel. 
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While the strong phase is vital to determining 7, the real part of the amplitude is still useful in and of itself. This 
is obtained from direct computation on the lattice from four-point function studies or via approximation methods, 
directly yielding a lower bound on the absolute value of the amplitude, which would be a valuable constraint on 7 
extraction. Comparing with the progress on if — >■ 27r decays [22l. [23l| . it is reasonable to foresee a lattice calculation 
which may reach a precision on the real part of this amplitude on the order of 15-20% within the next five years. 
With the input of such lattice computations, and the progress in the analysis of the CLEO data for the charm decays 
[l^ . one can envisage the extraction of 7 with error around 10%. 

Regarding the approximation methods for lattice studies, we first examine the possibility of chiral expansion of the 
real parts of B~ — )■ D^P~ and B~ — >■ D P~ amplitudes, in the framework of heavy-meson chiral perturbation theory 
(HMXPT) which merges heavy quark effective theory and chiral perturbation theory (XPT) |24l - l29| . The presence 
of the b and c quarks, both heavier than Aqcd, allows for systematic expansion in terms of AqcD/m^.c ("i^^c is the 
6, c quark mass). This expansion has already been used both in lattice determinations to leptonic and semi-leptonic 
decays, as well as in the continuum. Combined with the chiral expansion, it leads to a powerful tool for extrapolating 
lattice data to the physical pion mass. This extrapolation will still be an essential step in lattice calculations in the 
foreseeable future, since most lattice simulations are not yet performed at the physical pion mass. 

To begin, we examine the validity of HMXPT for the processes we are interested in. In the limit where both the 6 
and c quarks are treated as static, resulting in a soft final-state Goldstonc boson, this approach is valid. However, 
this limit is far from the physical regime, and such an extrapolation would introduce significant systematic errors. 
Therefore, the straightforward applicability of XPT is questionable (i.e. it would be a poor approximation with rather 
large errors). On the other hand, if wc perform simulations near the physical kinematic point of the decays of interest 
[B — > DK{ti)], the emerging D meson and the Goldstonc boson arc hard, with p ^ 2 GcV."' 

The appearance of hard external momenta does not, as one may initially assume, lead to a breakdown in the chiral 
expansion. It was recently shown that treating the s-quark as heavy and using SU(2) XPT works quite well for chiral 
extrapolations (sol . |3]| . This method can be generalised to processes in which external pions have hard momenta. 



and applications have appeared in analyzing decays [32[, K — 27r |33j . as well as extensions to semi- leptonic 



_B-dccays [34| . A central concern in all of these applications is how well XPT works in the presence of hard momenta. 
In particular, there is evidence that the hard pion docs not spoil the chiral logarithms, at least to ncxt-to-lcading 
order (NLO) [H^ 

Applying this to processes involving D mesons is straightforward, and a key result is that the hard momenta of the 
external mesons (both the D and the pion) will be absorbed into a redefinition of the low-energy constants (LEC's), 
and thus all remaining quantities will be soft. Thus we can still treat the D meson using the non-relativistic approach 
of HMXPT, so that corrections arising in the D sector will arise at ©(Aqcd/A/c), as usual. 

To investigate the relevant B decay processes, we are interested in the following current-current, A6 = 1, operators 
(a, j3 are color indices) 

(5) 
(6) 

(7) 
(8) 

For the decay channels B^ — > D'-'K^ or B^ — > D"K^ , we will set ^ s and for B^ — >■ D'^n^ or B^ — )■ iD^n^ we 
have = d. The corresponding effective Hamiltonian for these decays is 

= % E E [K;K,.Q(/^)Qr^'' + K;v^c,.Q(A*)Qr''' + h-c-] • (9) 

V ^ j=l,2 i=d,s 



— 


(9^7^(1- 




75)"/j) , 




(9L7''(1- 


75)&/j)(c/37m(1 - 


lb)Ua) , 




(C7''(l- 


75)6a)(u/}7M(l - 


75)C/3) , 




(9^7^(1- 


75)6/j)(w/37A.(l - 


75)Cq) . 


K- or B- 




we will set = 


s and fo 



* Note that this impUes large discretization errors of the form (ap)", with n > 0, and as such it would require either very fine lattices or 

a choice of action which would largely suppress these errors. 
^ As noted in Ref. |33| . hard-pion xPT may not be applicable to the extraction of the imaginary parts of the K — > tttt amplitudes. We 

will comment on this issue for nonleptonic B decays in Sec. |V] 
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We will focus on the nonleptonic decays which have the underlying processes b — >■ cud, b — > cus, b — )- ucd, and 
b — > ucs. The first two will be mapped onto different operators in the chiral theory than the final two, because 
they belong to different irreducible representations under the chiral transformation. Furthermore, a chiral field which 
creates a heavy-light meson with a c quark is not a field which destroys a heavy-light meson with a c anti-quark. We 
will discuss the details of these operators in the chiral effective theory in Sec. IIIIl 

The outline for this paper is as follows. First, in Sec. |ll]we present an introduction to HMXPT for B,D and B,D 
mesons. In Sec. Illli we construct the XPT operators for the quark-level operators in Eqs. ([S])-®. We then treat 
the leading-order calculation of B — >■ DP and B — >• DP and relate them to the unphysical B ^ D and B ^ D 
processes. In Scc. lIVl we discuss the tree- level resonance (initially either a B* or D*) contributions to the nonleptonic 
B decays in the framework of HMXPT. Finally we sketch the steps required for a full one-loop calculation in Sec. IVl 
and conclude in Sec. |Vll 



II. HEAVY-MESON CHIRAL PERTURBATION THEORY 



The strong-interaction chiral Lagrangian for the Goldstonc bosons is (the rj' is already integrated out) [3^, |33| 

£g = yTr (a,,I]a^I]t) + l^/2Tr (A^S + MY.^) , (10) 

where /i is a low-energy constant (LEG) related to the chiral condensate, M is the light-quark mass matrix, 

M = A\&g{mu,md,'ms) , (11) 

S = exp(2i$//) is the non-linear Goldstone particle field, with $ being the matrix containing the standard Goldstone 
fields, and we use a normalization for / such that « 130.7 MeV. Under an SU(3)l ® SU(3)r chiral rotation, E 
transforms as 

E — > LY.R\ where L e SU(3)l, and R e SU(3)r. (12) 

To account for the light-quark dynamics in hea vy-hght mesons, one can combine the formulations for heavy quark 
effective theory (HQET) and XPT into HMXPT [24H29l|. There is a U{2m) spin-flavor symmetry on the heavy-quark 
side for m heavy quarks, and the standard (broken) SU(3)l ® SU(3)r, chiral symmetry for the light quarks. 

We will sketch the relevant details for constructing HMXPT, using the notation of Ref. [1^. First, we have the field 
which destroys (creates) a heavy-light meson 

- (^) (7-V;W) _ ^,p(Q)) , , /ff (Q)t^o ^ (^.^*(Q)t + ^^p(Q)t) (^i+i^ , (13) 

where a is the light quark flavor index, Q is the heavy-quark index, and v is the four- velocity of the heavy-light meson. 
We use V for the heavy-light pscudoscalar field and V* for the heavy-light vector field. For the heavy-light fields with 
heavy anti-quarks, we have [39| 

It is convenient when dealing with both charm and bottom quarks and antiquarks to combine them into multiplets 
which transform under the U{4) spin/flavor symmetry. 

Suppressing the light flavour and velocity indices, under the heavy-quark spin/flavour transformation 5 G U{4), and 
the unbroken light-flavour transformation U(a;), the above fields transform as 

HQ{x)^SHQix) Ut (x) , Hq (x) ^ l]{x) Hq (x) , 

H^ix) ^ U(x) H^{x) 5t , Hq{x) -> S Hq{x) (x) , (16) 



H. 
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where U(x) is a function of L, R and 

The Goldstone bosons couple to the heavy-hght mesons in HMXPT via the non-linear realisation 

a=ys = e'*/^ 

which transforms as 

a{x) -^La{x) U'' {x) = U(.t) ct(.t) i?^ , a\x) Ra\x) V'^ix) = U(a;) a'' (x) . 



(17) 



(18) 



Due to the properties of the heavy-light meson fields in Eq. p6p . it is convenient to define objects involving the a 
field that transform only with U and U^. The two possibilities with a single derivative are 



(19) 



(20) 



The Lorentz vector can be combined with the derivative to form a covariant derivative acting on the heavy-light 
field or its conjugate: 



He 



ha 



(21) 



with implicit sums over repeated indices, and similarly for the Hq ^ fields. The covariant derivatives and transform 
under the unbroken light-flavour symmetry as 



H^^ ^ (i?&^)Ut 



UA^Ut , 



where we have dropped all the indices for simplicity. 

The leading-order chiral Lagrangian is given by £lo = + 'Chl.i, where 



(22) 
(23) 

(24) 



Tr means the complete trace over light quark flavor indices, heavy quark flavor indices, and, where relevant, Dirac 
indices. Since H and H always appear together in the Lagrangian, we treat HH as a matrix in light-quark flavor 
space: {HH)ab = HaHi,. The axial coupling g^^^ in the above Lagrangian determines the Goldstone and 

D*—!)— Goldstone interaction strength. Its value, g^r ~ 0.45, has recently been computed using unquenched lattice 
QCD [SSI. 



At the next-to- leading order (NLO), the Lagrangian contains a number of additional terms [38|, |43|, |4J| . Among these 
terms, only one of them. 



AoTr 



f±-H, 



(25) 



is relevant to this paper (A2 is a LEG). This operator breaks the heavy-quark spin symmetry and results in the B*—B 
and D*—D mass splittings. Notice that M-p is taken to be the corresponding B and D meson masses in this work, 
and we do not include other effects related to the breaking of heavy-quark flavour symmetry. 



III. THE CHIRAL EXPANSION FOR B DK{tt) AMPLITUDES 



The difficulty in the use of XPT in computations for B — > DP decay amplitudes originates in the large momenta 
carried by the final state hadrons. In general, the chiral expansion is known to be applicable only to processes involving 
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momenta well below the chiral symmetry breaking scale. On the other hand, it has been established recently that XPT 
can be valid for amplitudes containing hard final state particles (32| - [3^ l45j . One important point in such procedures 
is that the LEC's are no longer universal quantities for a fixed number of sea quarks. Rather, they depend on the 
hard momentum scale which results from either the kinematics or the mass of the external particles. 

This procedure of separating the hard scales in a process is described in detail in the references given above. The key 
point in this separation lies with the derivative couplings that give rise to momentum dependence in XPT calculations. 
When these momenta are external and hard, they can be absorbed into the LEC's of the theory. We will discuss this 
procedure explicitly with an example diagram for the process B — > DP in Sec. [V] 

First we discuss the construction of the XPT weak operators corresponding to those in Eqs. ([5])- ([8]). Omitting the 
colour indices which do not play a role in XPT, these operators can be written as 



gf-^.' = {ql Fi b) {UL T2 c) , (26) 

where = d or s, and 

= (^) ^' 

Pi =F2 = ri= Fa =7^(1 - 75). (27) 

Under the SU(3)l ® SU(3)r chiral symmetry group, Q^^'^'^ is in the (8l, 1r) representation, while Q''^'^'^ is in the 
(6l, 1r) representation. To bosonise these operators, we promote ri^2 and Pi, 2 to be spurion fields which transform 



Pi ^ L Pi 5^ P2 ^ ^ F2 L\ 



f 1 ^ L Pi S\ f2 ^ i P2 S\ 



(28) 



under the heavy-quark spin/flavour and chiral rotations. This renders the operators in Eq. (j26p invariant with respect 
to such transformations. We then find the bosonisation results in the leading order (LO) operators 



^x.. = E {"i.^^TrD [(<Ti,i?!;;,) P2S;S,Pi (<,'4) 

X 

for Q''-"='\ and 



"2,0: TrD 



a2.2;Tri 



Trr 



2.ri ( 7?g4 



}, (29) 



(i?iVli)]TrD [s.Pi 



(30) 

for Q''^"'', where Pro means the trace in Dirac space, and the summation over repeated indices are assumed. The 
symbols Sf^, and S^, are all possible pairs of Dirac structures allowed by symmetries [46j , 



(31) 



{75, 75} , {7^75, 7''75} , {^'75, ^75} , 75} , {75, hs} 



In particular, the positions of these Dirac structures in HMXPT weak operators are constrained by heavy-quark 
spin/fiavour symmetry. They have to be inserted to account for light-quark and gluon dynamics. 
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(a) 



(6) 



FIG. 1: Tree-level diagrams contributing to (a) B ^ D and (b) B — s- DP at lowest order, with no insertions of the strong 
Lagrangian. The box is the weak operator, the solid line is a heavy- light pseudoscalar (either B or D), and the dashed line is 
the light meson P. 



Performing the Dirac traces in Eqs. and ([50]) . we obtain 



O 



(c)t 



t 



[/3i + (/3i+/32)(^''-t')] 
+ [(/3i-/32)t;'^-/3i«^] 



(32) 



where Pi is a linear combination of ai^x and a2^x while /3j is a linear combination of ai^x and q;2,x- At the lowest order 
in the chiral expansion, only the first terms in the above operators contribute to _B — > DP and B — > DP processes. 
It is straightforward to demonstrate that if we evaluate the diagrams in Fig. [l]at leading order. 



{D"K-\0^JB-) = {D\-\0^^d\B-) = j{D-\0^,s\B-), 



(^K-\0^,s\B-) = {d\-\0^,,\B-) - j{D-\0^,,\B- 



(33) 



From Eq. ([5^ . it is clear that beyond the LO, the chiral expansion may become very different for B D^P and 

B~ — > 1^ P~ amplitudes. In the next two sections, we will discuss the generic features of these amplitudes at the 
NLO and leave the details to a future publication. 



IV. RESONANCE CONTRIBUTIONS 



In this section we discuss one generic feature of _B — > DP correlators and amplitudes, namely, the resonance contri- 
bution.^ This is partly incorporated in HMXPT via the inclusion of the vector heavy-light mesons. Figure[2l^b) shows 
a typical diagram in which a resonance {D* in this case) appears in the B — )■ DP correlators. One can also include 
heavier resonances in the effective theory [13], but it is beyond the scope of this paper. Here we will address the issue 



^ The conclusion presented in this section is also valid for B — s> DP decays. 
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FIG. 2: Tree-level correlators contributing to B — DP. The box is the weak operator, with (a) being the direct B — DP term 
and (b) being the term with an intermediate resonance (here a D*). 



regarding the contribution from the resonance in the time- momentum representation of correlators. To avoid com- 
phcations arising from the formulation of HQET and HMXPT in Euclidean space [1^, we work in Minkowski space 
with the comment that we also carried out a similar calculations by modelling the heavy-light mesons as relativistic 
particles in Eulidcan space and obtained the same conclusions presented in this section. 

We first set up the calculation for the LO correlator in Fig. [D^a). To mimic the setting for most lattice calculations, 
we integrate over the spatial volume for the positions of external B, D and pion (kaon) sources/sinks, i.e., we perform 
a Fourier transform for the spatial directions for each of the external points. On the other hand, we fix the location 
of the weak operator (the square in the diagram) to be at the origin. To be consistent with the notation in Sec. IIIIl 
we denote the velocity of B and B* by v and that of D and D* by v' . For simplicity, the velocity v is chosen to be 



1,0 , 



(34) 



and the time-ordering is implemented as 



< < < <p. 



(35) 



where tsjjp is the temporal locations of the B, D, P mesons, respectively. Using the Feynman rules derived from the 
HMXPT Lagrangian and the weak operators in Eqs. and ([5^ . the result for the contribution from this diagram 
in the correlator is 



f 



'iujptp 



2ujp 



gsDP (l\ ( e-'*°*° \ /e-'"p^p 



(36) 



where 



5d = v'-pd, and LUp^JMj,+pj, 



(37) 



with pr) and pp denoting the spatial momenta of the D and the Goldstone boson. The coupling qbop is one of the 
linear combinations of the EEC's (3i in Eq. ([32]) . 



Next, we discuss the correlator depicted in Fig. EJb). This diagram is calculated by integrating over the entire 
space-time for the location of the strong vertex (denoted by the circle) . It leads to the result 




-0(tz3)e-*^"*" 



-iujptp 



2ujp 



^i(u)p+So-^np)tn _ \ 



2iVQ{ujp + 5d - ^DP) 



^i{ujp+Sn~^Dp)tn _ ^ 
2v'q{uJp + 'Sd - Azjp) 



(38) 



where 



Adp = v' ■ (pd +P7r) + —r^ 

^0 



(39) 
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N X 
» » 

n* 

FIG. 3: A contribution to the process DP — >■ DP involving a single resonance. 



with Ajj denoting the D*—D mass splitting resulting from the heavy-quark spin symmetry breaking term in Eq. ([25j). 
When the final-state momenta are tuned sueh that 

ujp + 6d = Adp, (40) 

the resonance is on-shell and the correlator contains a linear term in to. 




^^cs\^p^Sd=Adp ~ f [nil o„./ I ( o, ) ( o^./ ) ' (4^) 

which is an energy shift of the final state. When one takes the ratio of the B DP correlator and the square root 
of the DP DP correlator, the tjj dependence arising from the square brackets in Eq. psp (hence this linear term 
in tjj), is exactly cancelled by the contribution from the diagram in Fig. [31^ The coupling is defined in Eq. (|24l) . 
and gsD* is a linear combination of the EEC's /J, in Eq. (|32p . Notice that gsD' is different from gsop and thus 
the resonance contribution results in general in an additional unknown parameter for B — > DP amplitude at the tree 
level. 



V. BEYOND TREE LEVEL 



The use of tree-level XPT is certainly limiting. While many systematic errors should cancel when looking at the ratio 
of B — )■ DP to B — )■ DP, going to higher order both in the chiral and heavy quark expansions is essential. Using the 
symmetry relations resulting from Eq. p2p . we could attempt an NLO calculation to make similar relationships at 
higher order, which is possible in the case of iiT —> 27r, as in, for example, Refs. (s^ - IH^ ). 



In order to treat these processes in the physical regime, we use the methods of Refs. [32| - |34| |: Hard-pion XPT (HPXPT). 
As discussed earlier, this formalism uses the fact that one or more of the momenta in the final state very well may be 
hard, and at the physical point for B — >■ DP, this is true. For this section we will focus on P = tt. 

In order to apply HPXPT to both B — > Dtt and B — > Dtt, there are quite a few one-loop diagrams that we must 
evaluate. The result of a complete calculation [i.e., the sum of all one- loop diagrams) is expected to take the following 
generic form [working with the SU(2) chiral theory for now] 



1 



In 



A2 



(42) 



where M is one of the the particular amplitudes from Sec. IIIIl and A^''''^° its tree-level value, a is a coefficient 
that depends on the particular kinematics chosen for the diagram, and L is a linear combination of low-energy 
constants as well as terms arising from higher-order chiral-level weak operators.® These would be determined from 
evaluating the full one-loop corrections to these amplitudes. We stress that a and L above will depend on all the hard 



^ This cancellation may not occur in partially-quenched QCD due to the loss of unitarity [49l45lll . 

* Repeating the spurion analysis of Sec. [II] would show in principle roughly 3-4 times as many LEC's arising at NLO relative to LO, but 
only certain combinations arise in calculations, eg., Eg. 1421 and as such there will effectively only be a small number of LEC's. 
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FIG. 4: One of the many one-loop diagrams that contribute to B — >■ Dtt, specifically one which shows the essential features 
that arise in HPXPT. 



quantities, specifically the mass of the external D meson and the momenta of both the external D meson and pion. 
This dependence is not known analytically, and it makes the LEC's non-universal when varying the hard momenta. 
However, at any fixed kinematics, the values of the LEC's are still fixed^. Additionally, we note that since all the 
hard scales are absorbed into the LEC's, we expect similar convergence as that of ordinary XPT. Corrections to the 
heavy quark expansion will be more significant coming from the D-meson, and thus both a and L will have 0{1/Md) 
corrections. 

In order to understand the specific details, we work through an example diagram, shown in Fig. |4l To evaluate this 
diagram, we envision a lattice simulation where momentum will be conserved at the strong vertex, but need not be 
at the weak vertex. Thus, we define the momentum entering the weak vertex as Pwk, the momentum flowing through 
the pion line is £ (the integration variable), and the external D meson has velocity v' and residual momentum fc, so 
that this diagram has the form. 



{D°TT-\0^,d\B~y""' f d'^e i iv' ■ {£ -~ p^) 



-/ , (43) 

where the coefficient arises from the weak vertex, and the momentum injected into the weak vertex, Pwk, is related 
to those carried by the external B, D and pion, 

Pit+Pd^Pb+ Pwk ■ (44) 

A = Mo — 1\Ib is the D-B meson mass splitting (which is of order l/rric — \/mb) and fc^ is the residual momentum 
of the £?-meson. 



This integral can be evaluated simply to obtain 

v'-k + A 



167r2/2 



v'ik +p.) + A + "'•^'^ + + ^ + - 1^ [-A2 



(45) 



with 



/2(m,(5) = -2(52 In j^—j -4<52F(m/J) + 2,52 , ^^g^ 
, fVl -a;2tanh"Vl -a;2 , < .x < 1 

= \ m — -1 m — T ^ 1 ■ ^^^1 

— Va; — 1 tan \j x'^ — 1 , x>\ 



We can examine this case in the limit where v'-k ^ m-,^^ which is the hard-pion limit. In this limit we find 

n2 



/2(m^,z;'-(/c+p^) + A) « -m^ hi ( ^ j 



(48) 



In practical lattice calculations, one would have to vary the pion mass, and extrapolate to the physical point. In this procedure, it is 
inevitable to change the momenta, and therefore the values of the LEC's. Fortunately, since the hard momenta arc all much larger than 
the typical pion masses in present and future lattice simulations, changes in the latter will result in very small variations of the former. 
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so that the full integral contributing to this diagram becomes either 

2 / 2 \ 

^(P^-0)^-2Y^l^i(#j ' (49) 
if we insert momentum into the weak vertex such that w 0, or 

if we choose p^k such that p^^ k, k. These would give rise to different values of the coefficient a in Eq. (05]). This can 
thus be extended to all of the diagrams that would contribute to one-loop order, and for each chosen set of kinematics, 
wc would be able to find different expressions for a in Eq. (|42p . and in general, the LEG L in that equation would 
have an unknown dependence on the kinematics. However, the pion mass dependence is well determined. 

We close this section by noting that HPXPT is not applicable for extracting the strong phases of B decays 
via the computation of the one-loop diagram in Fig. 21 The imaginary part in this diagram is proportional to 
y/^^' ■ (Pti- + ~ '^wj therefore grows with the increasing momenta carried by the final-state mesons, leading to the 
failure of the chiral expansion when and p-^ are large. This can be understood by noting that the imaginary part 
arises from the contribution in which both mesons in the loop are on-shell, and therefore cannot be soft. 



VI. SUMMARY AND OUTLOOK 



In this paper, we proposed a strategy for calculating B — >■ DP and B — > DP {P is a Goldstone boson) decay 
amplitudes via lattice calculations. Indeed the real part is accessible directly via four-point function calculations on 
the lattice as it does not suffer from the Maiani- Testa No-Go Theorem, though it is computationally demanding. 
As an approximation, one can invoke the chiral expansion, specifically taking into account the large momenta of the 
final state mesons. We argue that this hard-pion chiral expansion is valid for these decays, for similar reasons to 
those in semileptonic B decays and in K ^ rnr amplitudes. In general, this hard-pion chiral expansion results in 
momentum dependence of low-energy constants and the coefficients of the chiral logarithms. From our investigation of 
the structure of a typical one- loop diagram (Fig.|3|), it is shown explicitly how this occurs for the B DP amplitudes. 

We constructed the leading-order operators, relevant to these decays, in the chiral effective theory. We studied the 
tree-level resonance contributions in the framework of HMXPT, and showed that these contributions are accompanied 
by combinations of the LEG's which are different from that for the corresponding leading-order B — > DP and B DP 
amplitudes. As such, incorporating resonances in the study of the lattice correlators allows us to extract some of the 
leg's that are not accessible by applying XPT naively. 

To complete this inital approach, the complete one- loop contributions must be calculated [55| . One can combine these 
XPT results with lattice simulations to compute the real parts of S — > DP and B — > DP decay amplitudes. Although 
the lattice calculation for the imaginary parts of these matrix elements is challenging, their real parts can already 
provide important information for an accurate determination of the angle 7 in the b—d unitarity triangle in the GKM 
matrix. The real part gives a lower bound to the absolute value of the amplitude, which would be very useful in the 
phenomenology of 7-extraction, and by combining this with the information on strong phases from B and D decays, 
the absolute magnitude of the amplitude can also be deduced. 

Finally, let us note that, in the long run, as the lattice program succeeds in evaluating r^p and with experimental 
studies using larger data samples, experiment will be able to pin down rsp with increasing precision. We envision 
that a combination of these efforts could lead to an improvement in determinations of 7 to about 10% in 3-5 years. 
In the longer term, with the use of even more powerful computers and with data from Super-LHGb and Super- i? 
factories, the error could be reduced to a few percent. These improved determinations should allow a useful constraint 
on Vq^^° [Eq. ([4])] and consequently on Vut since all the other factors therein are already known quite well. Given 
the serious difficulties j56| in a precise determination of Vut through the conventional semileptonic methods, having 
an independent constraint via purely hadronic decays: B — > D{D)P may also prove useful. 
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